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I .  Introduction 

Consider  a  (linear  or  nonlinear)  dynamical  system  with  a  wide-band  noise 
input.  It  is  often  of  considerable  interest  to  approximate  such  systems  by 
diffusion  models  so  that,  e.g.,  Markov  process  techniques  can  be  used.  In  Ill  - 
(4 1,  171,  several  powerful  methods  for  doinq  this  have  been  developed. 

Roughly,  the  input  noise  process  is  parametrized  by  i  and  as  t  ■+  0,  the  bandwidth 
(BW)  ■»,  while  the  power  per  unit  BW  converges  to  a  constant.  The  limit  process 
is  found  via  methods  of  weak  convergence  theory.  The  methods  are  particularly 
useful  when  the  system  noise  (and/or  signal)  is  processed  nonlinearly;  i.e., 
only  nonlinear  functions  of  the  noise  appear  in  the  dynamics.  The  problem  is 
often  not  what  the  so-called  correction  term  might  be,  but  what  the  entire 
form  of  the  limit  is,  and  this  is  not  usually  easy.  In  fact,  when  nonlinear 

functions  of  the  noise  appear,  the  notion  of  "correction  term"  loses  much  of  its  sense. 
In  this  paper,  the  system  of  Figure  1  is  dealt  with. 

vE  »  F(vE)  +  DyE, 

u.i> 

y  ”  le  sign  u,  u  =  s  +  n  -  G(v),  v  (t)  6E  R  , 

where  nC(*)  is  a  scalar-valued  wide-band  noise  input  process.  Conditions 
on  F(  • ) ,  ne(*)  and  G(0  will  be  given  below.  The  main  result  is  that  as  e  -*■  0 
(BW  -*■  <■) ,  the  measures  of  { ( • ) }  converge  to  those  of  v(-)  where  v(*)  satisfies 
the  ltd  equation 

(1.2)  dv  -  F(v) dt  +  LD[  (S~~-— ? )  (^2A)  dt  +  h  In  2/a  dB], 

where  B(’)  is  a  standard  Brownian  motion  and  L,  a,  a  will  be  defined  below. 

Roughly,  "a"  is  related  to  the  correlation  function  of  the  ne ( •)/(E|ne(t) | 2) 

2  r 

and  o  is  the  intensity  of  the  spectrum  of  n  (■)  in  any  band  (0,BW)  for  siaall  c. 
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If,  in  the  system  of  Fig.  1,  the  saturator  and  gain  L£  were  replaced  by  a 
gain  K,  then  the  limit  would  be  (follows  from  the  method  of  this  paper) 

(1.3)  dv  «  F(v) dt  +  KD [ (s-G (v) ) dt  +  /2R  db], 

where  for  the  noise  model  used  for  (l*i)  ,  r  *  o  /a.  There  is  no  so-called 
"correction  term".  Owing  to  the  form  of  the  saturator  function,  the  formal 
technique  of  Stratonovich  is  inapplicable. 

The  example  is  offered  to  illustrate  what  can  be  done  with  one  particularly 
annoying  but  useful  nonlinearity.  The  basic  method  is  widely  applicable.  The 
scheme  is  unrelated  to  statistical  linearization,  which  in  fact  is  not  con¬ 
cerned  directly  with  approximating  processes. 

Before  proceeding.,  compare  (1.2)  and  (1.3)  for  the  case  when  the  feedback 
-G ( • )  is  supposed  to  be  stabilizing  (i.e.,  when  the  system  is  designed  to  make^-^^ 
the  error  s (t) -G(vE (t) )  small.  In  (1.2),  the  term  in  the  dynamics  whichr^ 
involves  the  error  is  proportional  to  l/o,  and  in  (1.3),  the  noise  term  is 
proportional  to  a.  Thus  for  small  o,  we  expect  the  limiter  to  enhance  stability 

without  increasing  the  noise  effects,  an  important  point  to  note.  For  large  o, 
the  limiter  does  not  seem  to  be  helpful.  A  simulation  comparison  of  the  "pre-limit" 
with  the  limit  for  a  somewhat  different  problem  (a  phase-locked  loop  with  a  satur¬ 
ator)  suggests  that  the  limit  (c  0)  results  are  often  the  "worst"  case,  in  that 
(for  example)  the  limit  mean  square  error  often  increased  to  the  limit  value  as 
c  -*•  0.  (They  also  suggest  that  often  the  limit  process  is  approached  quite  fast 
(as  measured,  say,  by  the  mean  square  value  of  the  input  to  the  limiter)  as  f  -*•  0.) 

Me  do  not  know  the  extent  of  applicability  of  this  rule  -  but  it  seems  to  hold 
frequently.  When  it  does  hold,  the  limit  results  can  provide  useful  upper  bounds, 
and  system  improvements  suggested  by  the  form  of  the  limit  might  well  be  improvements 
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for  the  “pre-limit"  case  also.  Unfortunately,  it  is  not  usually  possible  to  get 
approximate  diffusion  processes  where  the  BW  is  not  large  -  so,  even  if  it  is  not 
large,  the  results  for  large  BW  might  bo  a  useful  guide  to  the  qualitative  behavior. 

Reference  15]  contains  some  applications  to  problems  in  communications 
theory  of  the  same  general  idea.  But  owinq  to  the  unbounded  nature  of  the  noise 
and  the  form  of  the  discontinuity  and  feedback,  the  problem  here  is  harder 
and  the  analytical  details  different. 

Section  11  gives  specific  assumptions.  The  main  background  theorem  and 
some  comnents  on  weak  convergence  appear  in  Section  III,  and  the  convergence 
{ve(  •)  }  -*■  v ( • )  is  proved  in  Section  IV.  A  'similar  method  would  be  used  with 
other  nonlinearities. 


II.  Model  Assumptions 

1.  The  noise  model.  Let  z(*)  denote  a  stationary  Gaussian  process  with 
correlation  function  o2  exp  -a|t|,  a  >  0,  and  set  nE(t)  ■  zE(t)/e,  where 

z  (t)  «  z(t/e  ).  As  e  -*  0,  the  spectral  density  of  n  (•)  converges  to  2a  /a 
on  any  finite  interval.  The  scaling  is  a  convenient  and  common  way  of  getting 
a  noise  process  nE(")  whose  spectrum  converges  (as  e  •+  0)  to  that  of  a  white 
noise  with  a  constant  power/unit  bandwidth.  For  other  correlation  functions 
the  In  2  in  (1.2)  is  replaced  by  something  slightly  different.  We  use 
the  noise  form  only  to  facilitate  the  evaluation  of  the  coefficient  of  dB(*) 
in  (1.2).  The  Gaussian  assumption  simplifies  the  proof  that  certain  integrals 
converge  -  but  is  not  essential. 

2.  The  limiter  gain  Lc.  If  t£  ‘  L,  a  nunber  not  depending  on  e,  then 

as  c  ►  0,  the  "increased  wildness"  of  nE(*)  essentially  wipes  out  the  saturator  - 
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rep lacing  it  by  an  open  circuit.  Thus  L£  must  increase  as  c  decreases.  In 
any  particular  fixed  practical  system,  one  particular  value  of  L£  will  be  used. 
But  as  the  bandwidth  co»  this  value  of  L£  will  have  to  increase  (see  proof  in 
Section  IV)  and  eL£  will  have  to  converge  to  a  non-zero  number.  So  we  use 
-  L/e. 

3.  Other  assumptions.  G(-),  F(0  are  continuously  differentiable  and 
the  solution  to  (1.2)  is  unique  in  the  weak  sense.  s(‘)  is  right  continuous 
and  uniformly  bounded  on  [0,*) .  The  method  is  most  easy  to  use  if  the  functions 
are  smooth .  The  analysis  will  be  done  with  90 ( • )  replacing  g ( * )  »  sat ( • ) , 
where  the  piecewise  linear  90(")  is  defined  in  Fig.  2.  We  then  get  the  result 
{ve(-)}  -*■  v(-)  as  e  -*•  0,  then  o  0. 


III.  Weak  Convergence;  A  Convergence  Theorem 

X  X 

Tightness .  Let  D  [0,«)  denote  the  space  of  R  -valued  functions  on  [0,®) 

which  are  right  continuous  and  have  left-hand  limits.  A  certain  topology 

called  the  Skorokhod  topology  ((6],  section  14)  is  usually  put  on  Dr.  The 

process  v  (*)  is  considered  to  be  a  random  variable  with  values  in  Dr[0,<»)  and 

induces  a  measure  P  on  it.  {P  }  or  {vc(*)}  is  said  to  be  tight  iff  for 
e  e  —  *  ■■ 

each  6  >  0  there  is  a  compact  K. £  Dr[0,»)  such  that  P  (K.)  >  1-6,  all  e. 

o  £  0 

£  y 

{v  (•)}  is  said  to  converge  weakly  to  a  process  v(’)  with  paths  in  D  {0 ,«•)  and 
inducing  measure  P  on  it  iff  for  each  bounded  real-valued  function  g(’)  on 
Dr[0,<»),  / q(w)dP£(w)  /q(u)dp(u)  as  e  -*■  0.  Thus  weak  convergence  is  a 

generalization  of  convergence  in  distribution,  it  is  the  appropriate  form  of 
convergence  for  our  problm.  The  tightness  condition  for  (v£  ( •) )  will  hold 


under  our  assumptions 


where  t>N<v)  >  1  for  v  £  »  {v:  |v|<N},  bN<v)  =  0  for  v  £  SN+1  and  bN<v)  £  [0,1] 

and  has  third  derivatives  that  are  bounded  uniformly  in  v  and  N.  If  we  can 
prove  convergence  for  {v  '  (•),  z-*0)  for  each  N,  then  Theorem  1  says  that  we 


can  prove  it  for  (1.1).  Thus,  the  truncation  is  purely  technical  and  does  not 


affect  the  result. 


Definitions .  Let  A  denote  the  infinitesimal  operator  of  the  diffusion  (1.2). 
Let^'**  denote  the  a-algebra  induced  by  (vE'N(s)  ,nEfs)  ,s<t)  and  ee,N  the  cor¬ 
responding  conditional  expectation.  Actually  j?£'N  and  EE,N  depend  on  a  also. 

But  we  usually  suppress  the  a  affix.  Let  &  be  the  class  of  measurable  (u,t) 
functions  such  that  if  g(*)G  HI,  then  E|g(t+6)-g(t)  |  ■>  0  as  6  +  0  and 
sup  E|g(t)  l  <  «  and  g(t)  depends  only  on  (vE ,N(s) ,nE (s) ,s<t) .  We  say  p-lim^j  f^  » 
0  iff  sup*.  .Ejffi(t)  |  and  E|f^ (t)  |  0  as  5  0.  Define  an  operator  AE,N  and  its 

Oft 

domain  as  followss  g€  .<3?(AE'N)  and  AE,Ng  -  q  iff  g,q E  t#  and 

EE'Ng(t+fi)-g(t) 

p-lira  E| - - - q(t)  |  -  0. 

6|0  6 

The  following  theorem  is  Theorem  1  of  [2],  adapted  to  our  case.  denotes 
the  set  of  continuous  real-valued  functions  on  Rr  *  (0,») . 
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Theorem  1.  Let  the  equation  (1.2)  have  a  unique  weak-sense  solution.  Fix  W. 
For  each  f ( • )  £  Q)  ,  a  dense  set  (sup  norm)  in  <§q ,  let  there  be  a  sequence 
(fe'N(-) } £  &  satisfying  the  following; 


(3.2)  p-lim  |fc,N(t)-f(v",N(t)  ,t)  |  ->  0, 

e-K) 

tv*-0 

(3.3)  fC'V)  €  £^(At'N)  , 


(3.4)  p-lim  |AE,Nfe'N(t)-(AN-»|^)f(ve'N(t)  ,t)  |  =  0, 


e-K) 

a-*0 


.N  . 


where  A  is  the  infinitesimal  operator  of  some  diffusion  process  and  the  coeffic- 
snl 

{ ve  ( - ) }  v(-)  weakly . 


ients  of  AN  and  A  are  equal  for  v  £  SN_  Then  if  {vE'N(-) }  is  tight  for  each  N, 


Comment .  Tightness  is  not  hard  to  prove  here.  See  comments  at  the  end 

of  the  proof  of  Theorem  2,  which  applies  Theorem  1  to  our  case  (1.1).  Given 

f(*,*),  the  main  problem  is  to  find  the  fe,N(’)  and  to  verify  (3.2)  -  (3.4) 

(and  ultimately  to  prove  tightness).  The  method  used  here  and  in  [1],  [2]  is 

similar  to  the  averaging  method  used  in  (3].  We  choose  the  form  f  '  (t)  ■ 

f(ve,N(t),t)  +  fE,N(t)  +  fE,N(t),  where  fE'N(t)  is  chosen  so  as  to  "average  out" 

certain  noise-dependent  terms  in  Ae'Nf <ve,I,(t) rt) ,  and  fE,N(t)  is  chosen  to 

*C  N 

"average  out”  certain  noise-dependent  term  which  result  from  applying  A  ' 

£  N 

to  f , '  (t).  In  the  proof  lint  .  means  lim,^,  lim 

J.  €-*v  W  0*HJ 

a-K) 
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IV.  The  Convergence  Theorem 

Theorem  2.  Under  the  assumptions  in  Section  II,  {vt(-)}  converges  weakly  to 
v  ( • )  as  e  ■>  0  and  then  a  ->  0 . 

Proof .  Let  ^  ^  1 ,  the  subspace  of  fe’0  of  functions  whose  mixed  partial 

derivatives  up  to  order  2  in  t  and  3  in  x  are  continuous.  By  Theorem  1,  for 

each  N  and  f  (•,•)€  we  only  need  to  find  { f C ' N  <  - )  }  satisfying  (3.2)  -  (3.4). 

£  N  C 

For  notational  convenience  write  v  *  ( ■ )  as  v  ( •)  in  this  proof,  but  we  are  always 

e  N 

working  with  the  truncated  process  v  '  ( • ) , 

Part  1.  Fix  fC,*)  £  Then 

(4.1)  A€,Nf (vE(t) ,t)  =  ft<vE(t),t)  + 

+  bN(vE(t))f^(vE(t),t)  •  [F(vE(t)>  +  —  ga(s(t)+nE(t)-G(vE(t)>)). 
Note  that  for  u  in  any  bounded  set 

(4.2)  ^£g  (u+n£(t) )  s  gE(u)/c  =  ~[P{z(0) >-eu+ea}-P{*<0) <-eu-ea} +  0(ea)/e] 

g  a  a  e 

=  /2/ir  %  +  0(a)  +  0(c) , 

which  justifies  the  L£  *  L/e  scaling.  We  will  get  fE,N(*)  in  the  form 
fC,N(t)  -  f(vE(t),t)  ♦  f[,N(t)  +  f j '^(t) , 
where  the  fE,N(*)  will  be  defined  below. 
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The  following  estimate  will  be  used. 


On  the  set  {jz(t)|<l)  or  even  on  { |  z(t)  |^eaT//2} , 

-a.T 

|p{z(t+T)  €  B  |  z( t) }  -  p{ z ( t+T )  (  B}|  <  Ce  for  some  constants  C  and 
a^  >  0,  uniformly  in  B.  Similarly,  on  the  same  z(t)  set  and  for 


T2  >  T1  *  °' 


"Vl 


|p{z(t+t.)6  B.,  i—1 , 2 1  z(t)  }  -  P{z(ti+t)F  B.  ,  i*l,2)|<Ce  for  some 


a^  >  0  and  C  <  •  and  all  B^ ,  B2< 


In  the  sequel  the  values  of  a^  and  C  may  change  from  usage  to  usage. 

Define  g  (u,n£(t))  ■  g  (u+nE(t))  -  Eg  (u+nE(t))  and  define 

a  a 


fE,N(t)  -  -  b(vE(tn  [f  (vE(t)  ,t+T)DEE,Ng  (s(t+T)-G(veU>)  ,nE(t+t))dT 
1  e  N  J  v  t  a 

0 

ao 

•  eU>N(vc  (t) )  Jf  Mve(t)  ,t+e2t)DEE  ,Nga(s(t+E2T)-G(vE(t)  ),z(^+T)/e)dT 


By  (4.3),  fE,N(t)  “  0(e)  uniformly  in  w,  u,on  the  set  { | z(t/e2) | <1 ) . 
Define  w^  -  min{i:  e  aT/,2|  z(t/e2)  |<1} .  Write  fE'N(*)  as 


fE,N(t)  -  e 


w  « 

J1  Ee'V)  dT  +  ej  Ee,N 


{ • )  dt . 


The  first  ten  is  bounded  in  absolute  value  by  cCw^  and  the  integrand  of  the 
second  by  C  exp  -a^T.  Thus 


(4.4)  |fJ'N<t) |  <  Cc(l+w1)  <  Cell  +  max(0, log |zE(t) |) ] 
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Part  2.  It  can  be  verified  that  fE,N(0  €■  ^(Ae,N)  and  that 


Lb  ( v  ( t) ) 

(4.5)  A£'Nf^'N(t)  - - - -  fMvt(t),t)Dgn(s(t)-G(vC(t))<ne(t)) 


+  -  ]E^'NIb  (vC  (t) )  f 1  (v*  (t) ,  t+x)Dg  (s(t+i)-G(vl  <t)) ,nC(t+t)) ] ’  vL(t), 

6ft  N  V  (X 

o 

where  the  subscript* v  denotes  the  gradient  of  the  bracketed  expression  with 

respect  to  v€(t).  At  this  point,  let  us  simplify  the  notation  by  dropping  the 

b.  (v)  terms.  All  of  the  f.'  will  be  proportional  to  either  b.„(v  (t))  or 
N  -  1  N 

b2(vE(t)).  Changing  variables  T/e2  t  and  splitting  the  integral  in  (4.5) 

N 

6  £  £ 

into  two  parts  and  using  v  =  bN<v  ) [s-G(v  )+Dga/e]  (but  dropping  the  bN(v)) 
yields 

oo 

(4.6)  L2jEE'N(D'f  <vE(t)  ,t+e2T))  *Dg  (s(t+E2T)-G(vt'(t))  • 

0 

ga(s(t)-G(vc(t))+z(t/e2)/e)dT  +  o(e) 

OO 

+  L2fEE'N(D'fv(ve(t)  ,t+€2T))g{;v(s(t+e2T)-G(vf  (t)  ,z(~+x)/e)D  • 

0 

ga(s(t)-G(vE(t))+z(t/e2)/e)dT  +  0(e). 


7tie  terms  in  (4.6)  exist  by  the  same  arguments  which  led  to  (4.4).  He 
next  show  that  the  second  integral  of  (4.6)  is  negligible  as  e  0,  o  •+  0  and 
get  an  estimate  which  is  useful  for  the  tightness  argument.  The  facts  that 
a(t)  and  vc(t)  are  bounded  (recall  that  we  are  using  the  truncated  process 
O.D)  and  that  the  support  of  g„  (u)  is*  in  [-a, a]  and  that  |g  (u)  I  <  C/a 

*gauC- )  is  the  derivative  of  gaf*I  with  respect  to  its  argument.  The  subscript  v 
denotes  the  derivative  with  respect  to  the  explicit  argument  vs  replace  vC  (t^'take 
the  derivative  with  respect  to  v,  then  set  v  ■  vc(t). 
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will  be  used  frequently  and  perhaps  without  specific  mention.  Let  I { A }  denote 
the  indicator  function  of  the  set  A. 

By  (4.3)  it  can  be  verified  that 


(4.7)  Y  =  |E^'Ng  (s-C(v)  ,z(-~ +T)/r.)  |  <  [exp  -a  t  +  I {  |  z ( t/c^| >eaT//2 )  )C/e  . 

t  Cl  |  V  1 


We  need  a  bound  on  Y  which  goes  to  zero  as  e  >0.  First  we  get  such  a  bound 
when  |z(t/e2)|  _>  1.  Note  that 


(4.8)  P{ |s-G(v)+z(^-+t)/e j  £a  z(t/e2)=z  }  =  O(ue) 


uniformly  for  | zQ |  ^  1  and  T  >  0  (recall  that  s,  v  are  in  a  bounded  set  - 

for  each  N) .  Now,  (4.8)  and  the  facts  cited  above  (4.7)  imply  that  Y  is  bounded 

by  0(e),  uniformly  in  |z(t/e  )|  >.  1,  T  >_  0.  Thus  on  |z(t/e  )  |  >.  1, 


(4.9)  Y  <  [exp  -axT  +  X { | z( t/e 2) | >eaT/2> ] 1/2C (e/a) 1/2 

(use  | x |  <  a,  |x|  <_  b  |x|  £  /ab) .  Thus,  on  integrating  the  bound  when 
2 

|z(t/e  ) |  >_  1,  we  see  that  the  second  term  of  (4.6)  is  bounded  above  by 

C [1  +  max(0,  log|z(t/e2)  |)  ]  (e/ot)1//2. 

2 

Now,  we  look  for  a  bound  when  |z(t/e  ) |  <  1 .  Split  the  second  integral  in 

e  °° 

(4.6)  into  the  two  parts  /  +  /.  The  first  part  is  0(e/o) .  Note  that  the 

t  0  e  2 

density  of  z(-^+t),  t  t,  conditioned  on  any  value  of  |z(t/e  )|  in  [0,11,  is 

c  ,/2 

bounded  above  by  0(l//e) .  So  (4.8)  then  holds  with  0(ca)  replaced  by  0(ex/  a). 


'll' 


Conbining  this  estimate  with  (4.7)  yields  that  Y  is  bounded  above  by  (4.9)  when 
|z(t/e2)|  <_  1,  but  with  the  change  that  le/a)1^2  is  replaced  by  (c1//4/ a1//2)  in 

(4.9) .  Thus,  on  integrating  the  bound,  we  get  that  the  second  term  of  (4.6)  is 
bounded  above  by 

(4.10)  C 11  +  max(0,  log |  z( t/t2)  j )  ]  . 


Part  3.  We  turn  our  attention  to  the  first  term  of  (4.6)  and  show  that, 
by  an  "averaging",  it  can  effectively  be  replaced  by  its  expectation.  To  facil¬ 
itate  the  development,  we  define  the  following  terms. 


h  (v,t,x,p)  =  L2D'f  (v,t+T+p)D-g  (s(t+T+p)-G(v)+nE(t+T+p) )  • 
e  w  a 


•  g  (s(t+p)-G(v)+n  (t+p)), 
01 


I*  /  .  \  -  T  2  >  _  /  .  2  .  2.  %  ^  _2_.  2_  v  „ 

H  (v,t,T,p)  =  L  D  f  (V, 


(4.11) 


(4.12) 


...t+e  t+e  p)D-g  (s+(t te  t+c  p)-G(v)+z(— r+T+p)/e) 
w  a  i 

c 


•  g  (s(t+e2p)-G(v)+z(-^+p)/e)  , 
a  £2 

CO  00 

A£,Nf(v,t)  H  |eh^ (v,t,T,0)dx  -  ~  |EhE(v,t,T,0)dr. 

0  E  0 
oo  oo 

f2'N(t)  S  ^2  |dp  j  dT[E£'NhE<vE(t),t,T,p>  "  =hc(v,t,T,p) 

£  0  0 
00  00 

jdp  |  dTlE£'NHE(vE(t) ,t,T,p)  -  EHe (v, t, T ,p) 


v»ve(t) 


0  0 


w  (t) 


], 
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where  I  implies  that  v  is  replaced  by  vF(t)  after  taking  the  expectation. 

'v  ve(t) 

C  N 

We  must  show  that  (•)  is  well  defined.  First  note  that  the  inner  integral 
£  N 

of  f2'  (t)  (with  p  =  0  and  a  change  of  variables)  is  just  the  first  term  of 

r  N 

(4.6)  centered  about  its  expectation.  The  form  of  f  (•)  is  chosen  to  allow 

us  to  average  out  the  first  term  of  (4.0)  and  to  effectively  replace  it  by  its 
f  N  £ 

average  value  AQ*  f(v  (t),t). 

By  the  method  used  to  bound  |fF'N(t) |,  we  can  get  that  the  inner  integral 
of  fF'N(t)  exists  for  each  p.  Recall  the  definition  =  min{w:  e-aw/2 | z( t/s2) | <1 } 
and  write  (4.12)  as 

OO  03  tt> 

C2  I  dp  j  dT  EF,NB  +  e2  |  dp  j  dr  EF,NB  =  II  +  I. 

0  0  0 

First  we  show  that  II  is  well  defined.  By  (4.3)  and  the  definition  of  w^,  the 
absolute  value  of  the  integrand  in  II  is  bounded  above  by  C  exp  -a^p,  a^  >  0. 

Also  |eH£ (v, t,T ,p) |  C  exp  -a^r  for  some  a^  >  0.  By  (4.3)  and  on  the  set 
{p>wx) ,  and  for  C,  a^  (whose  values  again  may  change  from  usage  to  usage) 

|EF'NHe(v,t,T,p) |  <  EF,N|eE'N2  H  <v,t,T,p> | 

t+E  p 

<  C  EE'N[exp  -a^T  +  l{e’aT/,2|  z(-^+p)  |^l)l 

*  C  exp  -ajT  +  CP{  |z(-^+p)  |>eaT//2j  pHr^zU /e2) } 

e 

<  C  exp  -a^  ♦  Ce’aT//2E{  |z(^j+p)  I J  p>Wj »z(t/e2)  } 

<  c  exp  -a^i . 
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Chebychev's  inequality  is  used  to  get  the  next-to-last  inequality.  Conbining 

the  above  estimates  yields  that  the  integrand  in  II  is  bounded  by  (for  sene  a^  >  0, 

2 

C  <  <»)  C  oxp  -a1<t+p).  Hence  II  “  D(r  ). 

2  2 

The  term  I  is  also  0(e  )  but  not  uniformly  in  z(tA  ) •  Bound  the  inner 
integral  of  I  by 


OO  00 

J  d,|E‘'1V|  <  ECt-"  |  dT  |Ee’N 


t+e2p 


8  III. 


By  the  arguments  used  to  get  the  bound  on  |fE,N(t) | ,  we  get 


III  <  Ce2  EE,N[1  +  max(0,  log|  z(-~+p)  | )  ] 


<  Cc2  Ee,N[1  +  log(|z(-^+p)  |  +  1)J. 


By  Jensen's  inequality  and  the  concavity  of  log(’). 


III  <  Ce2  (1  +  log(  |  z(^r-)  |  +  C)  ] 


2 

Since  Wj  <  C  max(0,  log|z(t/e  )|), 

(4.13)  | f|#N{t> |  <  Ce2[l  +  log(|*(ij) |  +C>]2. 

E 

Henceforth,  we  will  give  only  an  outline  of  the  details,  which  can  all 
be  filled  in  via  the  estimates  and  techniques  developed  above.  It  can  be  shown 
that  fE,N(*)  €  #(Ae,N)  and  that 
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(4.14)  *e,Nfe,N^tj  =  negative  of  first  term  of  (4.6)  +  A^'^f (vC(t) ,t) 


♦  (terns  whose  i>-lim  equal  zero), 
i  K) 

(t  KD 

The  term  whose  p-lim  =  0  is  just  (flj’^ft))  *vL(t),  where  fG'^  *s  5*'*<ii®nt 
of  the  expression  for  fG,N  with  respect  to  the  argument  v*  (t) .  The  components 
of  Ae'NfG,N  which  involve  are  bounded  by  0(e).  Loosely  speaking,  the 

remaining  component  is  of  the  form 


(4.15)  o(e)  +  E|{tEt,Nfw^a,v9a  +  Et'Nfw9a9«,v  “  E  fw9a,v9a 


-  E  f  g  g  ]drdp  *  g  . 
w  a,v  a 


where  we  omit  the  function  arguments.  By  a  method  similar  to  that  used  to  get 

2 

(4.13),  we  get  the  bound  (4.13)  on  (4.15)  but  with  (e/a)  replacing  e  . 


Part  4.  The  estimates  obtained  in  Parts  1-3  imply  that 


p-lim  |fe(t)-f (ve(t) ,t) |  *0, 
e-*0 
o-K) 

p-lim  |Ae'NfG(t)-ft(v€(t),t)  -L/2A  A«(t)tSlv-(t)2j  D'fv(vC(t) ,t) 


AG,Nf(ve(t),t) |  «  0. 


A  proof  very  similar  to  that  in  ([5],  Section  6,  part  2)  yield/  that  AG'**f(v,t) 
D’f  (v,t)D(ln  2) /a  uniformly  in  v  for  each  t.  In  calculating  the  limit,  the 


*0ne  of  the  seasons  for  the  choice  cov[*(0) ,s(t)]  -  o2axp  -at  is  to  allow  us  to  save  j 

mock  by  using  this  result,  lhe  choice  allowed  an  explicit  evaluation  of  the  diffusion 
tmm.  With  other  choices  the  diffusion  coefficient  would  be  left  in  an  "integral"  form,  j 
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G(  • )  ,  b ( ' )  play  no  role  and  the  limit  (e  ■*  0,  a  ■+  0)  is  the  same  as  for  the 

case  (a  -  0,  e  -*•  0)  .  If  the  b  (v)  terms  were  retained,  the  result  would  be  the 

N 

2 

same,  except  that  either  b  or  b  .  would  multiply  the  f  ,  f  or  f  .By  what 

N  N  v  w  wv  1 

has  just  been  said 

(4.16)  p-lim  |Ae,Nfe,N(t)  -  <|r  +  A)f(vE(t) ,t) |  -0, 

e-*0  3t 

a-»0 

where  A  is  the  infinitesimal  operator  of  v(*)  in  (1.2).  If  b  were  retained, 

N 

the  A  in  (4.16)  would  be  replaced  by  some  AN  which  would  equal  A  where  b  (•)  *1, 
i.e.  in  SN>  Thus,  by  Theorem  1,  if  (vE'N(-)}  were  tight,  then  the  proof  would 
be  completed. 

Tightness .  Use  (12],  Theorem  2). The  conditions  of  Theorem  2  [2]  hold  if  (4.17) 
holds  for  each  N  and  T  <  «•; 

(4.17)  lim  Tim  P(sup|AE'NfE'N(t) |  >  K>  -  0, 

K-k»  e-*0 

a-*0 

lim  P{sup  | fE,N(t)+f2*N(t)  |  5}  ■  0,  each  6  >  0. 

e-K)  t<T  1  i 

o-K) 

But  (4.17)  follows  from  (4.10),  (4.13)  (and  a  similar  estimate  for  (4.15),  and 
the  fact  that  the  Gaussianness  and  stationarity  imply  that  («<-  T  >  c>^ 

lim  sup  eY|z(t/c2)|  »  0  w.p.  1.  Q.E.D. 

e-K)  t<T 
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